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Non-commutative geometries have been proposed as an approach to quantum gravity and have 
led to the construction of non-commutative black holes, whose interior singularities are purportedly 
eliminated due to quantum effects. Here we find evidence that these black holes are in fact unstable, 
with infalling matter near the Cauchy (inner) horizon being subject to an infinite blueshift of the type 
that has been repeatedly demonstrated for the Reissner-Nordstrom black hole. This instability is 
present even when an ultraviolet cutoff (induced by anticipated non-commutative geometric effects) 
to a field propagating in that spacetime is included. We demonstrate this by following an analogous 
argument made for Reissner-Nordstrom black holes, and conclude that stability is dependent on the 
surface gravities k_ and k+ of the inner and outer horizons respectively. In general if k_ > k+ , as 
we show to be the case here, then the stability of the Cauchy horizon becomes highly questionable, 
contrary to recent claims. 



The development of a convincing and self-contained theory of quantum gravity is arguably the most important and 
difficult problem in theoretical physics today. While no such theory has yet achieved this feat, a general tenet thought 
to be characteristic of quantum gravity is the presence of a natural small-scale cutoff of spacetime; a "graininess" if you 
will. This cutoff is expected to arise when approaching the Planckian scale. One of the primary applications expected 
of quantum gravity is a better understanding of curvature singularities predicted by classical general relativity; such 
divergent behavior is indicative of a breakdown in the classical theory and it is thought that such extremes will be 
placated by a proper treatment using quantum gravity. As such, much work has been performed towards understanding 
black holes in the context of quantum gravity theories. Several metric solutions have been derived using arguments 
from these theories that are free of the r = singularities that plague classical solutions [IH5] . Indeed, when one 
imposes a small-scale cutoff one should hardly expect to find a divergence which classically results from trying to 
describe an infinitely small point in spacetime. 

Despite being regular at r — however, these quantum gravity black holes generically have more than one horizon. 
One therefore ought to be wary of their overall stability since it is well known that classical black holes with more 
than one horizon also exhibit divergent behavior at their inner horizon (the Cauchy horizon). This can be seen to 
arise from the infinite blueshift of external radiation that occurs on the horizon. The example of prime interest here 
will be the Reissner-Nordstrom black hole [6HH] , though this instability is also found in Kerr and Kerr-Newmann black 
holes. In this paper we investigate whether or not this instability is present in the black hole solution inspired by 
noncommutative geometry [3J \5\ . 

One proposed solution to quantum gravity, motivated by the principle of small-scale cutoff, is to impose a noncom- 
mutative structure to the spacetime manifold 9| . The mathematics of noncommutative geometry is a generalization of 
differential geometry and its application to spacetime results in an uncertainty relation of the form As^As" > 
where Q^ v is a matrix that quantifies the magnitude of the noncommutative deviation from a classical spacetime 
manifold; its entries are postulated to be on the order of the Planck area. This uncertainty relation implies that in 
this formulation there is no notion of a point-like event in spacetime, in the same way that in quantum mechanics 
there is no notion of an exact point in phase space. In addition to entirely changing the notion of what a manifold 
is, such a small-scale cutoff also changes the form of quantum fields propagating on the manifold by imposing a 
natural momentum cutoff on the order of 1/v 9, where 9 is defined to be the mean value of the matrix entries. Both 
of these modifications are important in predicting the phenomena arising from this theory. As with other quantum 
gravity theories, work has progressed in deriving black hole solutions inspired by noncommutative geometry [H[S]. 
The solution corresponding to the classical Schwarzschild metric is 

ds 2 = -f{r)dt 2 + fir^dr 2 + r 2 dn 2 
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The mass function m(r) is proportional to the total mass M and a lower, incomplete gamma function 7(0, b) = 
Jq x a ~ 1 e~ x dx. In the classical limit 6 — > we see that m(r) — > M and Eq. d 1 b reverts to the Schwarzschild metric. 
This metric admits two horizons provided M > 1.9-\/6\ below which point it becomes extremal. An example of this 
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is shown in Fig. (JlJ where we plot the metric function f(r) with 8 = 1 and the two cases M = 1.9 and M = 4. The 
limiting value of approximately 1.9 is determined from the conditions /(ro) = f'(vo) = for some value ro (the radius 
of the extremal horizon). The condition /'(ro) = can be used to find ro; one can show that it is given by ro = 
where x solves the equation 7(3/2, x) = 2x z / 2 e~ x . The value x is numerically determined to be x « 2.284, with which 
the other condition /(ro) = can be used to obtain Mj\J~Q = ^pKx~ x e x / '4 ~ 1.904. We will only consider conditions 
under which two horizons are present, and we will label the outer and inner horizons r + and r_ respectively. 





r 

(a) M = 1.9-v/0 (b) M = 4^0 

FIG. 1: Metric function /(r) 

A preliminary investigation into the stability of this black hole claimed that it was stable with respect to blucshiftcd 
radiation (with an ultraviolet cutoff) at the Cauchy horizon, since the energy density of an ingoing held pulse (specif- 
ically, a Gaussian pulse since a ^-function pulse is forbidden by the ultraviolet cutoff) was found to remain finite upon 
reaching the Cauchy horizon |10) . However this investigation neglected to examine perturbative scatterings from the 
primary pulse that occurs near the horizons. We will show here that such scattered waves to first order do have 
diverging energy densities at the Cauchy horizon, a result which, while not conclusive, is highly indicative of overall 
instability. Such a calculation was performed in [6j to indicate the instability of the Reissner-Nordstrom black hole, 
and our analysis here will follow from the same line of reasoning. 

We will find it convenient to use the tortoise coordinate r* and null coordinates u and v, defined as 



dr* = /(r) 1 dr 

t, v = r* +t 



(2) 
(3) 



The non-commutative structure of spacetime suggests that there is an ultraviolet cutoff to matter fields, and 
consistent with previous work [10] we shall impose the same on a massless matter field <fr by including a factor of 



2 e/4 



in the mode expansion. Neglecting angular dependence, this modified field takes the form 

00 duo 



where tp u solves 



and the potential is given by 



dr* 2 



V(r*(r)) 



9/A e -iut 
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(uj 2 -V(r*))^ 

f(r) f i(£ + l) 
r V t 



f(r) 



(4) 



(5) 



(6) 



where £ represents the fields angular mode. V(r*) vanishes at the horizons and near them can be shown to decay as 



V(r*) cx e 2K + r * 
V(r* 



for r — > r i or r* — > — oo 



-2«_r* 



for r — > r_ or r* — > +oo 



(7) 
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where k± are the surface gravities at the two horizons, defined as 

«± = ±i/'( r ±) ( g ) 

We are principally interested in the field near the inner horizon r_ since this is where stability comes into question. 
ipuj will in general consist of two linearly independent solutions corresponding to left-moving (ingoing) and right- 
moving (outgoing) waves traveling along surfaces of constant v and u respectively. Thus, at the inner horizon the 
field will be of the form 

— \g { -\u) + g {+ \v)\ (9) 
r_ L J 

We will need to derive the form of </ - ' and g( + ^ eventually. Assuming that they are given, however, we can compute 
the energy density p of the field as measured by a freely falling observer near the horizon with four- velocity U a ; we 
have p = <j>, a (t>,f)U a U p + \<i>,a4>' a ■ However, since u, v — const are null surfaces, the form of <j> near the horizons implies 
that this will be dominated by the |</> jQ [/ Q | 2 term. 

The four-velocity of a timelike, radial geodesic is found to be 

U * = J> U r = -^/E^f (10) 

where we define U r to be negative since this will always be the case between the horizons r_ < r < r + (since r 
necessarily decreases for any observer in this region). E is a constant of the motion, the sign of which gives the 
direction of travel between the horizons; E > corresponds to a left-moving observer and E < to a right-moving 
one. 

Our goal is to compute the energy density p oc \U a g[a %> \ 2 - It is easily seen that — ±gW and g,r = f~ 1 g^ ', 
where g^ is the derivative of g^ with respect to v or u depending on the sign. With this we obtain 

V a g<£> = £?-{±E-\I? - ft'*) (11) 

Let us examine this result in the limit r — > r_. Recalling that / vanishes at the horizons, for a left-moving 
observer (E > 0) we see that U a g^ remains finite but U a g\ a ' diverges, unless of course g(~) can compensate for the 
divergence. Conversely, for a right-moving observer (E < 0) U a g^ diverges while U a g] a remains finite. Assuming 
the black hole is non-extremal the metric function near r_ is / ~ — 2K_(r — r_), and so the divergence goes as 
U a g^ oc g^> (r — r_) _1 for r 

Ik— (r— r_ ) ' 



From the form of f near r_ we also obtain the tortoise coordinate: dr* ~ ^ — ^ — r, thus 



r* ~ -(2K_)" i ln|r-r_| (12) 

For a right-moving observer near r_ it is easily shown that = jjt ~ 1 , from which we obtain 

— v = —t — r* ~ — 2r* + const ~ kZ 1 log \r — r_ | + const 
=^> (r - r-)- 1 ot e K ~ v (13) 

This gives us the form of the divergence in U a g^ expressed in null coordinates (recall that as r —t r_ , v — > oo for 
a right-moving observer and u — > oo for a left-moving one). That is, 

U a g^ oc g (+) 'e K - v for r ~ r_ and E < (14) 

Thus, if the inner horizon is to remain stable then </ + ' must decay at least as fast as e~ K - v in order to placate 
the divergence as v — > oo. A similar analysis shows that g^ must decay at least as fast as e~ K - u in order to stop 
the divergence of U a g[ a ' as u — > oo for observers with > 0. Our goal now is now to compute these quantities to 
determine stability. 

In order to compute g^ we reproduce a calculation used in [5] to determine the inner horizon stability of the 
Reissner-Nordstrom black hole. Since we are only interested in the field near the horizons, where the potential is 
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exceedingly small, we can decompose the total solution to Eq. (pi), which we will call <J/ W , into the zero-potential 
solution plus an infinitesimal perturbation produced by the small potential: ^ u = ip u + e^. We set the 
initial conditions to consist only of ingoing waves and we choose the time dependence to be e~ lult ; this requires that 
tp u = e~ tur . An ingoing (left-moving) wave will scatter off of the small potential near the horizons, and these 
scatterings are represented by e w . In the analysis performed in [10] only the the main waves ipui were considered, and 
it was found that the energy density due to these waves remains finite on r_. However they neglected to consider 
scattered waves e w , and it is these waves that (as we will show) give rise to instability, for reasons similar to that for 
a Reissner-Nordstrom black hole [5J. 

There are two scatterings that are of potential interest with regard to stability at the inner horizon. Fig. [2] shows a 
section of the Penrose diagram (including the external universe and the r_ < r < r + region) with the two scatterings. 
The first of these consists of right-moving waves traveling along the left branch of r_ ; these would have scattered off 
of the main wave as it neared r_ and they are labeled with a 1 in Fig. For these waves we must check the form of 
g(' to test for stability. The second scattering of interest consists of left-moving waves which travel along the right 
branch of r_. These waves can form by the following process: consider a scattering produced just after the main wave 
has entered from the outer horizon; it would be right- moving and traveling along r + . As this wave approaches the 
intersection of r + and r_ in the Penrose diagram it enters a region of strong potential. Seeing as the wave magnitude 
is assumed to be very small, a strong potential would be expected to scatter this wave again with effectively 100% 
efficiency; in this case the entire wave is scattered such that it is now a left-moving wave traveling along the right 
branch of r_ as labeled by a 2 in Fig. M. For these waves we must check the form of to test for stability. 




FIG. 2: Displaying the different scatterings off of the main wave which are of potential interest. 



In order to solve for 



e u we use a Green's function G^^r* ,y*), 




(15) 



Having chosen an e 



time dependence, we opt for the solution 




(16) 



Now, since ^ = + e w solves Eq. ^ while ip t 
that e u approximately solves 



solves the same equation with the potential set to zero, we have 




(17) 



the solution of which is 




(18) 
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We must pause now as the form of e u is not actually the held that we should use if we are considering a noncommu- 



tative spacetime. This is not surprising once one recognizes that the (^-function in Eq. ( 15 ) is not an allowed object 
given a noncommutative context. We can therefore not use the canonical Green's function (^^(r*, y*) as we have 
defined it here. We instead must use a modified Green's function G^e (r* , y* ) , presented in [T^], which is derived in 
the noncommutative framework. This modified Green's function takes the form 

G ue {r*,y*)=e e ^G M {r*,y*) (19) 

9 a 2 

The operator e s*-*^ acts to smear out point-like objects in order to accommodate the spacetime uncertainty of 
noncommutative geometry. When applied to a (5-function, for example, it outputs a Gaussian of width ~ 1/V&- 



Given the solution for G^oi 1 "* ,V*) in Eq. (16) we can rather easily compute the modified Green's function. For the 
r* > y* region we have 



G u e(r*,y*) = E 



1 ™ d 2n \ ( 1 



n\ dr* 2n J V 2iu> 



l_ e Mr'-v') y H_! ( 6u ?) 
h.i.i ^ <n\ v ' 



e iuj(r*-y*) 



,2\n 



2iu> 

and a similar result follows from the r* < y* region. In total we have 



2iui ^—^ n\ 

71 = 

e Mr*-V*) e -8" 2 (20) 



G^g(r\y*) = \ 2 \" . , . f. fl2 V (21) 



The noncommutative-modified perturbation field that we wish to use is then given by 

poo 

e ue {r*)= G ue (r*,y*)V(y*)My*W (22) 



Here we are interested in the scattering produced by the outer horizon potential V(y*) = Voe 2K+v , y* —> — oo. For 
purposes of computational ease let us set V(y*) — for y* > 0. We evaluate this integral for r* > since this will 
be the case when the wave approaches r_, which is what we are interested in. Identifying ij} u = e" iLjr we obtain 
trivially 

e ue (r*) = I 06 , . - (23) 
Including time dependence we obtain the right-moving wave: 

t/ -iujUp — Ouj 2 

e-^eir*) = . (24) 

4lU)(K-{- — lUJ) 

We now need to evaluate the total wave consisting of all modes. Given e~ lut ipu = e~ %uv the main wave consists of 
an ingoing, Gaussian pulse (j> = J ^-e~ u e ' i e~ wt ip u — ^ry/je~ v ^ '■ The total perturbation over all modes is 

ee = - e-^e-^ed.^ 1 - V °f. ' e . du, (25) 
r J^^ r J_ OQ - iuj) 

As explained above, when this wave approaches the intersection of r + and r_ in the Penrose diagram it will be 
entering a region of high potential, which can be expected to scatter the small wave with near 100% efficiency. We 
assume this condition here, and so the wave after this second scattering will be of the same form as that just derived 
but traveling leftwards instead of rightwards. Recalling Eq. ([9]), we thus have the form of g^ + \v) for this wave 

V ° A e , £ s du, (26) 



6 



In order to test for stability, however, we require the derivative 



9 (+Y (v) = ^ ~ 



did 



(27) 



which integrates to 



(28) 



56* (k 4 



2v s 
59, 



2, so it is unimportant in the decay 



As v — > oo the complimentary error function approaches erfc 

behavior. Therefore, we see that (v) decays as e~ K+v for an observer approaching the right branch of r_. Recall 
that the necessary condition for stability is that (v) decays at least as fast as e~ K ~ v , and thus the problem now 
becomes one of comparing the two surface gravities n±. That is, if K- > k + then (v) does not decay fast enough 
to suppress the divergence present in Eq. ( 14 ) and r_ is unstable. For the noncommutative black hole this is indeed 
the case, as we will show. 

First recall that there were two possible divergences that could occur at r_ , the other one generated by fields near 
the left branch of r_ and with stability contingent on g^ (u) decaying at least as fast as e~ K - u . In this case a similar 
analysis can be performed and it can be shown that it does decay fast enough, and so the Cauchy horizon is stable in 
this respect. 

Turning to the more interesting case, the surface gravities are defined to be k± = ±|/'(r±), where the metric 
function f(r) is given by Eq. 0. Taking the derivative of / and using the horizon relation 2m(r±) = r± we find 



K± = (1 - 2m'(r±)) 

2r ± 

The metric function m(r) is given by a lower, incomplete gamma function, as in Eq. fl). 
easily computed, thus giving us expressions for k± as functions of r± 



(29) 



As such, its derivative is 



m (r) = 



M 



K± = ± 



20F6» 3 / 2 
1 



2r± 



1 



-r 2 /49 



M 



(30) 
(31) 



The task now is to compute the values of r± . Looking at Fig. ( lb ) (plotting the metric function for 9 = 1 and 
M = 4) we see that k_ > k + (that is, the slope is steeper at r_ than it is at r+), yielding an instability at r_ . To 
confirm this result in the general we consider two regimes: M » y/d and M w 1.9\/#. In the intermediate regime 
between these values we find that k_ > k + , either numerically or simply by plotting f(r), as in Fig. (lb). 



We first examine the M ^> s/9 regime. To do so we note that r± are given by the intersections of the mass function 
m(r) and the line r/2; to visualize this we plot the two functions for the case of 9 — 1 and M — 8 in Fig. pjl). We 
see that in this case r_ is small and r + is equal to 2M to a very good approximation (since the mass function levels 
off to the constant m(r) ~ M by the time r/2 intersects it). As we increase M with respect to \fd these statements 
only become stronger. 

Thus, in this regime we will take r + = 2M for the outer horizon. For the inner horizon we utilize the fact that 
r_ is small and we Taylor expand m(r) about r — 0; it can be shown that the first non-vanishing derivative is 



d m(r) 



dr 3 



r=0 



and so for small r the mass function is well approximated by a cubic function 



m(r) 



M 



6Vtt9 3 



r 3 + 0(r 5 ) 



(32) 



We compute r_ by finding the intersection of this with r/2; along with our result for r + we have 



2M and r_ 



M 



for M > Ve 



(33) 
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FIG. 3: A plot of the mass function m(r) and the line r/2 for the case of 9 — 1 and M = 8. The intersections 

determine the values of r±. 



Inserting these results into the surface gravities, Eq. (31), we obtain 

3 



> + 




1 - 



K - ~ 2^ 3yW2 ( 




(34) 
(35) 



from which it can be seen that k + will generally decrease and k_ increase as Mj\[Q increases. Indeed as M/yff — > oo 
we have k + —> 1 /AM and k_ — > oo, as should be expected. Having observed that K- > k + for intermediate values of 
M/Ve we thus conclude that this inequality will continue to hold as M/y6 increases indefinitely. 

Finally, we wish to examine the M 1.9v^ regime; in this case we simply evaluated r± numerically for several 
values of M and 9 and inserted these values into Eq. (31 1 to compute k±. The results of this exercise are displayed 
in Fig. Q where we plot k± versus Mj\f9 for 9 — 0.25 and 9 = 1. As can be seen, k_ continues to be greater than 
k+ as M/V9 -> 1.9. 




(a) 



0.25 




(b) 6 = 1 



FIG. 4: k± versus M/V9 for two values of 9. The upper and lower curves correspond to k_ and k + respectively, 

thus showing that k_ > n + for M ~ 1.9y9. 



We thus conclude that k_ > k + in general, and this implies that the Cauchy horizon of the noncommutative 
geometry inspired black hole is always unstable with respect to the perturbative scatterings described. Since the 
scatterings are assumed to be small the fact that they diverge reflects a breakdown in our approximation, and so we 
can not yet say with absolute certainty that the Cauchy horizon is indeed unstable. However, such a divergence is 
certainly indicative of instability, casting doubt on previous claims for this black hole |10j . We emphasize that this 
type of analysis is easily performed for any two-horizon black hole, although the scalar field modifications (if any) will 
be particular to the theory in question. For example, the authors have also applied similar reasoning to a black hole 
motivated by loop quantum gravity . The results of this work suggest that the loop black hole may have stronger 
stability properties than the one being scrutinized in this letter. 

One might argue that the infalling pulse of scalar field is not a generic form of infalling matter. We therefore 
consider a stress-energy corresponding to an ingoing null dust present on the manifold, but we neglect the back 
reaction this induces in the metric. This is of course an approximation, but no more of one than was already being 
assumed in the precious analysis. We will see that the resulting energy density as seen by a timelike observer crossing 
the Cauchy horizon will, instead of being of the form p ~ e 2K-v e -2n+v ^ j-^^gj. g like p ~ e 2K - v L(v), where L(v) is 
the "luminosity" of the null dust. L(v) is expected to follow an inverse power law in the large v limit, as first shown 
by Price |13j . Thus in this scenario there is no way to avoid diverging energy density as v — ¥ oo, unlike in the previous 
analysis where the competition between two exponentials afforded some hope of this. 

First we cast the metric, Eq. ([T]), in (v,r,6,<f>) coordinates; the result is 

ds 2 = -f(r)dv 2 + 2dvdr + r 2 dtt 2 (36) 

where the metric function f(r) is unchanged. 

We then consider a stress-energy tensor which corresponds to an ingoing null dust. It takes the form 

T aP = n(r,v)(d a v)(d fi v) (37) 

With these last two equations we determine the form of p,(r, v) by computing the conservation condition = T Q ^. 
The result is that we must have p(r,v) — A(r)L(v), where L(v) is an arbitrary function of v and A{r) must satisfy 
the equation 

= r dA ^+2A(r) (38) 
dr 



the solution of which is 



A(r) cx \ (39) 
Thus, without loss of generality (since L(v) is arbitrary) we conclude that p,(r,v) is given by 

tir,v) = 1 ± s L{v) (40) 

Far from the black hole the spacetime is well approximated by a Schwarzschild solution, and we note that this form 
of /u,(r, v) is exactly the same result as that obtained for the Schwarzschild black hole. We thus identify L(v) with the 
usual luminosity, and we expect the analysis by Price |13j to hold to a good approximation in our case as well. That 
is, we have L(v) ~ v 1 ~ p with p > 12. 

In order to determine the energy density a timelike observer would measure from this dust we also need to specify 
the four- velocity of the observer. This was previously done for TP and U r , and it is easy to show that U v takes the 
form 

TP = J^(E- VE 2 fir)) (41) 

where we recall that E < for an outgoing observer between the horizons; this is the type of observer that we are 
interested in. 

The energy density measure by such an observer is given by 

p = T a pU a UP=n(r,v)(U v ) 2 (42) 
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As the Cauchy horizon is approached we have f(r) — »• 0, and so the negativity of E implies that (U v ) 2 ~ AE 2 / f(r) 2 as 
this approach takes place. We also know that near the Cauchy horizon the metric function goes like f(r) ~ (r — r_ ) ~ 
e~ K - v , and so we obtain that the energy density has the behavior 

p ~ e tll - v L(v) (43) 

near the Cauchy horizon. Since follows an inverse power law, we conclude that the observer measures diverging 
energy density upon crossing the Cauchy horizon. 

Of course the preceding null dust calculation is not without its limitations. Apart from the neglect of back-reaction 
effects, it is much less obvious how to include noncommutative modifications to this form of infalling matter, unlike 
the situation for the scalar field pulse. However the latter is not a generic form of infalling matter. 

Notwithstanding such considerations, our calculations provide strong evidence that the noncommutative black hole, 
while retaining finite curvature at r — 0, may have an unstable Cauchy horizon. A more complete treatment of this 
problem will entail an analysis of the phenomenon of mass inflation [§] , in which the mass function of the Reissner- 
Nordstrom black hole diverges at the Cauchy horizon. This approach takes into account back-reaction effects from 
the geometry. The suggestion that mass inflation cannot occur for the noncommutative black hole with ultraviolet 
restricted fields [TU] appears to be predicated on the claim of finite energy densities, a result that we have now seen 
does not ensure stability. A more rigorous analysis of mass inflation in this case is needed before one can make 
conclusions either way. 
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